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Abstract — An optical transmitter that uses entangled light 
generated by spontaneous parametric downconversion (SPDC), 
in conjunction with an optimal quantum-optical receiver (whose 
implementation is not yet known) is in principle capable of 
obtaining up to a 6 dB gain in the error-probability exponent 
over the optimum-reception un-entangled coherent-state lidar to 
detect the presence of a far-away target subject to entanglement- 
breaking loss and noise in the free-space link [1], [2]. We present 
an explicit design of a structured quantum-illumination receiver, 
which in conjunction with the SPDC transmitter is shown to 
achieve up to a 3 dB error-exponent advantage over the classical 
sensor. Apart from being fairly feasible for a proof-of-principle 
demonstration, this is to our knowledge the first structured design 
of a quantum-optical sensor for target detection that outperforms 
the comparable best classical lidar sensor appreciably in a low- 
brightness, lossy and noisy operating regime. 

I. Introduction 

An optical transmitter is employed to interrogate a distant 
region engulfed in bright thermal light, suspected of containing 
a weakly reflecting target. The return light is processed by a 
receiver to decide whether or not the target is present. Recently 
[1], Lloyd, building up on work by Sacchi [4] showed that 
in the above scenario, a "quantum illumination" transmitter, 
i.e., one that uses entangled light at the transmitter and an 
optimal quantum receiver, can perform substantially better 
than an un-entangled coherent laser transmitter, despite there 
being no entanglement left between the target-return and the 
idler beams due to high loss and noise. This is the first 
example of an entanglement-based performance gain in the 
bosonic -channel setting where the initial entanglement does 
not survive the loss and noise in the system. More recently 
[2], Tan et. al. showed that using a sequence of identical two- 
mode-squeezed Gaussian states obtained from spontaneous 
parametric downconversion (SPDC), in conjunction with an 
optimal receiver, one may obtain up to a factor of 4 (i.e., 6 
dB) gain in the error-probability exponent over an optimum- 
reception coherent-state radar in a highly lossy and noisy 
scenario. This optimal receiver can be abstractly expressed 
as a projective measurement that projects onto the positive 
eigenspace of the difference of the density operators of the 
multi-mode states of the target return and retained idler modes, 
under the two hypotheses - target absent, and H\\ target 
present. However, no structured receiver design is yet known 
that can harness any of the 6 dB performance gain. 



We present the design of a structured receiver, which despite 
being inferior to the (yet un-implemented) optimal joint- 
detection scheme, in conjunction with the SPDC transmitter is 
shown to achieve up to a factor of 2 (i.e., 3 dB) error-exponent 
advantage over the optimum-reception classical sensor in the 
high loss, low brightness, high noise regime. Our receiver 
attempts to directly measure the remnant phase-sensitive cross- 
correlations between the return-idler mode pairs, by mixing the 
return and the idler beams on a parametric amplifier and photo- 
detecting the output. Numerical evidence suggests that our 
receiver achieves the asymptotic error-exponent of the optimal 
separable measurement, hence indicating that any superior 
receiver would have to make a complex joint measurement 
over multiple return-idler mode pairs. 

II. Background 
An SPDC transmitter generates K independent spatio- 
temporal signal-idler mode pairs |a^,6^|; k € 
{1,...,K}, each prepared in an identical entangled two- 
mode-squeezed state with a Fock-basis representation 
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where Ns is the mean photon number in each signal and idler 
mode. In the quadrature representation, \4>)si is a maximally- 
entangled zero-mean Gaussian state with second-order quadra- 
ture field moments given by 
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standard quadrature field commutators [ojy , afj] = and 

[asla { s}} = [4i!> a J» 1 = iWmn applyfor 1 < k < 
Each signal mode is transmitted serially toward a spatial region 
that may or may not contain a weakly-reflecting specular 
point target, but in either case, contains a bright thermal-noise 
bath. Under hypothesis Hq (no target), the target-return mode 

'We drop the superscript (fe) from the annihilation operators whenever 
convenient. As each signal-idler mode pair is prepared identically and undergo 
an identical channel transformation, this will cause no ambiguity in notation. 



Q>r — Q-B, where clb is in a thermal state with mean photon 
number Nb 1. Under hypothesis Hi (target present), cir — 
y/nas + \/l — KO.B, where the overall channel transmissivity 
K <C 1, and clb is now in a thermal state with mean photon 
number Nb/(1 — k). Under Hq the joint return-idler state for 
each of the K mode pairs p^j is a product of two zero-mean 
thermal states (pn b <S> Pn s ) with mean photon numbers Nb 
and Ns respectively. Under Hi, the joint return-idler state for 
each of the K mode pairs p^J is a zero-mean Gaussian state 
with second-order quadrature field moments given by 
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The binary detection problem is to decide between hy- 
potheses Hq and Hi (assuming they are equally likely) with 
minimum probability of error, using a quantum joint-detection 
measurement on the state of the 2K return-idler mode pairs at 
the receiver, (p^j)^ K ^ ', where h = or 1, given hypotheses 
H or Hi respectively. 



Helstrom derived the minimum probability of error P 
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for discriminating two quantum states as a function of the 
number of available copies (or trials) K [5], P 
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yl — J2 n 7» + ' j /2» where 7^ + ' ) are the non-negative eigenval- 
ues of - (Pri) { ® K) - Building up on the classical 
Chernoff bound in classical detection theory [6], recently 
Audenaert et. al. derived the quantum Chernoff bound (QCB) 
as an upper bound to P^ in , and also showed the QCB 
to yield the exact asymptotic rate exponent of the min- 
imum error probability [7]. Combining that with a (rela- 
tively loose) lower bound on P^ in [2], and defining Q s = 

Tr ((/ 5 i?/) s (/ 5 fl/) 1 " s ) and Qqcb = min < s <i Q s , we have 
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where the second inequality (QCB) is asymptotically tight 
as K — > oo. The QCB is customarily represented as 
Pemin — b~ KRq /2 in terms of an error-rate exponent Rq = 
— ln(QqcB)- The third inequality is a loose upper bound 
known as the Bhattacharyya bound. 

Symplectic decomposition of Gaussian-state covariance ma- 
trices was used to compute the QCB explicitly [3], for both 
the coherent-state and the entangled (SPDC) transmitters [2], 
and it was shown that in the high loss, weak transmission 
and bright background regime, i.e., with Ns <C 1, K <C 
1, and Nb ^> 1, the entangled transmitter yields a QCB 
error-exponent Rq = kNs/Nb, which is four times (or 6 
dB) higher than the error-exponent Rc — kNs/(4:Nb) for 
a coherent-state transmitter with identical per-mode average 
transmitted photons as the entangled transmitter. In Fig. [T] we 
plot the regions captured by the upper and lower bounds in |7]) 



GO 

- 





it) & 



M 
o 




-1 
-2 

-3 
-4 
-5 
-6 
-7 
-8 
-9 
-10. 



Coherent-state 
transmitter bounds 



Quantum-illumination 
transmitter bounds 




Classical transmitter & 



Homodyne Receiver 



Quantum illumination 
transmitter & OPA Receiver 



K = 0.01 

N s = 0.01 
N B = 20 



logio(*0 



Fig. 1 . Regions between quantum Chernoff upper bounds (solid curves) and 
(loose) lower bounds (dotted curves) on the minimum error probability jTJ, for 
coherent-state (shaded dark/blue) and quantum-illumination (shaded light/red) 
transmitters with K transmitted modes each with Ns = 0.01 mean transmit- 
ted photons per-mode, Ng = 20 mean thermal-noise photons per mode, and 
k = 0.01. The coherent- state transmitter lower bound also applies to ALL 
classical-state transmitters with a total of KN$ photons in all K transmitted 
modes; thus depicting an undisputed asymptotic advantage obtained by the 
quantum illumination (SPDC) transmitter with joint optimum reception over 
all K return-idler mode pairs. The curve plotted with (magenta) circles depicts 
the error-probability performance of the coherent- state transmitter and mode- 
by-mode homodyne-detection receiver, which is asymptotically optimal (as 
K — > oo) for the coherent-state transmitter in the Ng 2> 1 regime. The thick 
black line is the performance of the OPA receiver that we propose in this paper, 
which performs substantially better compared to the optimal classical sensor 
(magenta circles) — note that the magenta circles is the best error-probability 
performance that can be achieved by any classical transmitter (the error- 
probability lower bound of the blue region is rather loose and anything below 
the magenta circles cannot be achieved by any classical transmitter). The 
error-probability performance of the (yet unknown) optimum joint-detection 
receiver for the SPDC transmitter will lie in the (red/light shaded) region 
between the bounds for the quantum-illumination transmitter, and will hug 
the quantum-illumination transmitter upper bound for high values of K. 



for both the classical and the entangled transmitters, showing 
a clear advantage of quantum over classical illumination. 

III. Receiver design for classical illumination 

For a coherent-state transmitter, each received mode an 
is in a thermal state with mean photon number Nb, and a 
mean-field value (or) = or V kN$ for hypotheses Hq and 
Hi respectively. Hence, homodyne detection on each received 
mode Sjj yields a variance-(2A r B + l)/4 Gaussian-distributed 
random variable Xk with mean or t/kNs given the hy- 
pothesis. Assuming both hypotheses to be equally likelj|^] the 
minimum error probability rule to decide between H and Hi 
is to use the sufficient statistic X = Xi + . . . + Xk against 
a threshold detector, i.e., say "F " if X < {K^/k,N s )/2 

2 With no information about the prior probability of a target being present, 
we can safely assume the hypotheses to be equally likely. Unequal priors just 
change the inconsequential pre-factors in actual error-probability and Chernoff 
bound expressions presented in the paper. Error-exponents remain the same. 
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and "Hi" otherwise. The probability of error is given by the 
Gaussian error-function: 
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where the approximation holds for 
Rc bom = kN s /(AN b + 2) is the 
1, i? Chom ~ kN s /AN b , identical 
for the coherent-state transmitter, 
that for K large enough, P,, hom 
in the high-background regime, 
detection is asymptotically optimal 
mitter. 
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kN s K/(4N b +2) > 1 and 
error-exponent. For N B ~3> 
to the QCB error-exponent 
It is straightforward to see 
< e - KRc ^™/2. Therefore 
mode-by-mode homodyne 
for the coherent-state trans- 



IV. Receiver design for quantum illumination 

For the SPDC transmitter, each received return-idler 
mode pair {or, a/} is in a joint Gaussian state with zero 
mean under both hypotheses, and covariance matrix V = 

([ a R a/ a R a\ ] [ a R a\ clr aj ]) whose en- 
tries are readily computed using the quadrature field moments 
in Eqs. Q, <|3j and (|6j. For Hypothesis Hi, 
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and for hypothesis H , V = diag(A B + 1,N S + 1,N B ,N S ). 
Unlike in the coherent-state transmitter case, the entangled 
transmitter results in zero-mean joint return-idler states under 
both hypotheses. As is evident from the covariance matrices 
under H and Hi, in the Ns <C 1, k <C 1, N B ^S> 1 regime, 
the sole distinguishing factor between the two hypotheses 
that makes quantum illumination perform superior to the un- 
entangled coherent-state transmitter, are the off -diagonal terms 
of V bearing the remnant phase-sensitive cross-correlations 
of the return-idler mode pairs when the target is present, 
yj kNs(Ns + 1). The optimal joint-detection receiver acts on 
all the 2K received return-idler modes and yields the minimum 
probability of error P^ in by capturing the remnant return- 
idler phase-sensitive cross-correlations for the Hi -hypothesis 
in the most efficient way. 

The receiver approach we propose uses an optical paramet- 
ric amplifier (OPA) which is a type-II degenerate amplifier 
constructed from a non-linear crystal. The incident return 

(k) (k) 

and idler modes a R ' and a} ; are combined and amplified by 
an OPA driven by a strong local pump beam, producing pairs 
of output modes 
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where G > 1 is the gain of the OPA (see Fig. [2|. Each output 
mode c is in a zero-mean thermal state with mean photon 
number under the two hypotheses given by 

GN S + (G - 1)(1 + N B ) = N Q , H Q 



GN S + (G -l)(l + N B + kN s ) 
-2y/G(G- l)y/ K N s (N s + 1) 4 Ni, Hi 



{ 



TO 



Optical 
Parametric 
Amplifier 
(OPA) 

Gain, G > 1 



Total number of clicks registered 
in K received modes 



at(fc) 



Fig. 2. A quantum joint-detection receiver that mixes all the received K 
return-idler mode pairs pairwise on an optical parametric amplifier (OPA) 
with gain G, and counts the total number of clicks N on a photon counter at 
one output port of the OPA over all K output modes. The receiver decides 
in favor of hypotheses Ho or Hi depending upon whether N < iVth or 
N > iV t h, Nth being the decision threshold. 



i.e., p c = E I T=o(^™/( 1 + N m y+ n )\n){n\, for m G {0,1} 
for Hq and Hi respectively. Hence, the joint state of the 
K received modes &- k > is a A'-fold tensor product pf K of 
identical zero-mean thermal states with per-mode mean photon 
number N or Ni depending upon which of the two hypothesis 
is true. A A-fold product of thermal states is diagonal in 
the A'-fold tensor-product of photon-number bases of the K 
modes. Hence the optimum joint quantum measurement to 
distinguish between the two hypotheses is to count photons on 
each output mode and decide between the two hypotheses 
based on the total photon count N over all K detected modes, 
using a threshold detector. The probability mass function of 
Af under the two hypotheses is given by 
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where m = or 1 based on which hypothesis is true. 
The mean and variance of this distribution are KN m and 
Aof Tl respectively, where af n = N m (N m + 1). For large K 
the above conditional distributions for N approach Gaussian 
distributions (due to the central limit theorem (CLT)) with 
means and variances given by KN m and Ka^ respectively. 
The probability of error is given by 

1 



Pi.oPA = i^rfc (VPqpaA- 
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where an error-exponent Ropa = (Ni — No) 2 /2(<7o + fi) 2 
can be achieved using a threshold detector that decides in favor 
of hypotheses H or Hi depending upon whether N < N t ^ 
or N > N th , with N th = \K(ai N + a Q Ni)/(a Q + a{)] . 
Ropa is a function of the OPA gain G. Given Ns <C 1, 
K <C 1, Nb 1, intuitively a small gain, G = 1 + e 2 
with < e 1 will be optimal to distinguish between the 
two hypotheses, such that the difference in the total mean 
photon count is dominated by the term proportional to the 
phase-sensitive correlation, i.e., K(Ni — No) = K(e 2 nNs + 
2e V / rT^v / ^RiVsTT)) ps 2eKy/ K Ns(N s + 1), for 
Ns < 1, « < 1, N B > 1 and e < 1. For the problem 
parameters as chosen in Fig. [T] i.e., Ns — 0.01, k — 0.01 
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and Nb = 20, Popa is maximized for G = 1 + 5 x 1CP 3 , 
confirming our intuition that a small gain is optimal. For this 
gain, Popa = 2 x 10 -6 . For the above parameters, error- 
exponents for the classical and SPDC transmitters with optimal 
measurement are Rc — 1.25 x 10~ 6 and Rq — 5 x 10~ 6 
respectively. Hence for these parameters, our receiver gets as 2 
dB gain in the error-exponent over the classical system. 

Finally, to establish the asymptotic error-exponent per- 
formance of the OPA receiver, we will use the classical 
Bhattacharyya bound to the error-probability pj opA- The 
Bhattacharyya bound to the error probability in distinguishing 
between the two distributions P/v|# (rc|Po) and Pjv|ffi ( n \Hi) 
using K i.i.d. observations is given by: 
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with R B = {Ni - No) 2 /{SNo{l + No)). This translates to the 
bound 
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where the Bhattacharyya bound exponent Rb is given by 
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e 2 KNs(N S + 1) 
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kN s /2N b , (20) 



for a choice of e 2 = N s /VNb, for A s < 1, n < 1, 
A^ ^> 1 (e <C 1 is automatically satisfied J] As Popa > Pa, 
therefore by construction, for a weak transmitter operating in 
a highly lossy and noisy regime, the OPA receiver achieves 

3 A different e satisfying Ng/Ng <C e 2 <C l/Ng would work as well. 
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Fig. 3. Symbol-by-symbol (separable) Helstrom minimum error-probability 
measurement on each return-idler mode pair seems to have identical error- 
exponent as (and performs only slightly better than) the OPA receiver. 
The separable Helstrom measurement performance was calculated by ex- 
plicitly evaluating the single-shot minimum error-probability -P, ^in = 
^1 — ~^2 n 7^ + '^ /2, where 7,' + ^ are the non-negative eigenvalues of — 

p^\, using CLT to compute the error-probability with independent measure- 
ments on all return-idler mode pairs, followed by majority-vote hard decision. 



at least a 3 dB gain in error exponent over the optimum- 
receiver classical sensor whose QCB error exponent Rc = 
kN s /4N B - For N s < 1 and e < 1, both A and N ± < 1. 
Hence, a single-photon detector (as opposed to a full photon- 
counting measurement) suffices to achieve the performance of 
the receiver depicted in Fig. [2] 

V. Conclusions 

Entanglement has proven to be one of the most useful 
resources in quantum information, with applications to im- 
proving classical communication rates over noisy quantum 
channels, teleportation of unknown quantum states over long 
distances, and in quantum algorithms that can potentially 
solve certain problems (such as factoring) faster than the 
corresponding best known conventional classical algorithms. 
Whereas realizing useful quantum computing is not practically 
feasible with current technology, exploiting quantum effects 
to build optical communications and sensing systems that 
could perform better than corresponding classical systems 
seems to be well within the reach of current state-of-the-art in 
experimental quantum optics. 

"Quantum illumination" is a novel preliminary attempt at 
using an entangled source of light to detect the presence of 
a weakly-reflecting target in the far-field subject to a very 
lossy and noisy channel [1]. Even though the efficacy of 
using spontaneous parametric downconversion - a well-known 
entanglement source - as a transmitter to obtain a significant 
advantage (6 dB in error-exponent) over using a classical laser 
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transmitter for target detection was established by means of 
quantum Chernoff bounds [2], a structured receiver design that 
may harness this advantage is not known to date. Our work 
puts forth an explicit receiver design using conventional quan- 
tum optics, that when used along with an SPDC transmitter 
is capable of obtaining up to a 3 dB error-exponent gain over 
the optimum-reception classical sensor. 

There are several problems and challenges - both in theory 
and experiments - that remain open. It is still not clear whether 
a sequence of identically prepared two-mode-squeezed states 
constitutes the best quantum-illumination transmitter. It will 
hence be of interest to explore whether a transmitter that 
uses non-Gaussian entangled states, or complex higher-order 
entanglement between all the transmitted signal-idler mode 
pairs can achieve more than just a constant-factor improvement 
in the error exponent. It would also be interesting to extend 
this work to detecting and imaging spatially extended speckled 
targets. 

In terms of receiver design, the problem of constructing a 
receiver for the SPDC transmitter that achieves the full 6 dB 
error-exponent gain is certainly at the forefront. Numerical 
evidence suggests that any alternative receiver that makes 
independent measurements on each return-idler mode pair is 
not likely to perform significantly better than the OPA receiver 
in the absence of classical feedforward or soft-decision post- 
processing through the K measurement instance^] thus neces- 
sitating a superior receiver to make a joint measurement over 
all return-idler mode pairs (see Fig. pjr] The current model 
assumes the receiver to have complete knowledge of the signal 
power, channel loss and noise power for the optimal design 
of the receiver. Whereas the knowledge of the transmitted 
signal power at the receiver is a reasonable assumption, a 
full estimation-theoretic study of the performance of quantum 
illumination with no prior knowledge of channel loss and 
noise, would be necessary in order to design a realistic 
prototype of the quantum sensor. 
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Appendix A 

In order to substantiate our claim about the performance of 
single-shot Helstrom measurement (depicted in Fig. [3]), we 
compute the error-probability performance of the pair-wise 

4 The author thanks Baris I. Erkmen, JPL for pointing out that by phase- 
conjugating each received mode and detecting the conjugated return and 
retained idler modes on a balanced dual detector would achieve the same 
factor-of-2 improvement in the error-exponent as the OPA receiver. 

5 To see the details of the calculations, see Appendix A 



(separable) hard-decision Helstrom measurement (the quantum 
measurement that minimizes the single-shot probability of 
error) followed by a majority- vote detector after indepen- 
dently detecting all K return-idler mode pairs. In order to 
compute the Helstrom measurement minimum probability of 
error P^ in , we need to compute the non-negative eigenvalues 
of the difference of the density operators of the two-mode 
return-idler states under the two hypotheses, i.e. p^J — p^j. 
For our propagation and noise model, these two return-idler 
density operators can be computed in full generality in terms 
of their respective tensor-product Fock-state basis elements. 
For a pure-state two-mode entangled transmitter of the form 
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the return-idler state matrix-elements come out to be: 
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N b (N b + 1-k)_ 
mi — m = ni2 — ri2 is satisfied and 
i{m2\R{mi\p^ I \ni) R \n 2 ) 1 = for mi-ni ^ m 2 -n 2 . \n) R 
and \n)i for n € {0, 1, . . . , 00} are the photon-number states 
for the return and the idler mode respectively, and each form 
a complete orthonormal set of bases for the respective state 
spaces. 2F1 is the well-known hypergeometric function. Using 
the above expressions for the density operators, we computed 
the error-probability performance of the symbol-by-symbol 
(separable) Helstrom minimum error-probability measurement 
on each return-idler mode pair, and numerical results seem 
to suggest that this measurement has the same error-exponent 
as (and performs only slightly better than) the OPA receiver. 
The separable Helstrom measurement performance was calcu- 
lated by evaluating the single-shot minimum error-probability 
P e%n = f 1 - J2 n A where are the non-negative 



eigenvalues of p^j — p R ^, and then using the central limit 
theorem to compute the error-probability with independent 
measurements on all return-idler mode pairs, followed by a 
majority-vote hard decision (see Fig. [3]). 
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